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We calculate correlation function in the Einstein-Podolsky-Rosen type of experiment with massive 
relativistic Dirac particles in the framework of the quantum field theory formalism. We perform 
our calculations for states which are physically interesting and transforms covariantly under the full 
Lorentz group action, i.e. for pseudoscalar and vector state. 
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I. INTRODUCTION 



One of the most puzzling aspects of quantum me- 
chanics, its nonlocality, is illustrated by the Einstein- 
Podolsky-Rosen (EPR) paradox dj. Quantum mechan- 
ics predicts that for a pair of entangled particles, flying 
apart from each other, measurements give results incom- 
patible with our intuitive conceptions about reality and 
locality. Quantum mechanical predictions have been con- 
firmed in many EPR-type experiments. Most of these ex- 
periments have been performed with photon pairs. How- 
ever all EPR experiments with photons are subject to 
the so called detection loophole — low efficiency of pho- 
ton detection allows the possibility that the subensemble 
of detected events agrees with quantum mechanics even 
though the entire ensemble satisfies the requirements of 
local realism. Therefore the fair-sampling hypothesis, 
stating that the detected events fairly represent the en- 
tire ensemble, must be assumed. Detection loophole has 
been closed recently in the experiment with massive par- 
ticles [2J . This experiment in turn does not overcome the 
so called locality loophole in which the correlations of ap- 
parently separate events could result from unknown sub- 
luminal signals propagating between different particles. 
(Locality loophole has been closed in the experiment with 
photons Q.) One can hope that future experiments with 
relativistic massive particles could close both mentioned 
above loopholes. Therefore such experiments seem to be 
very interesting for the basics of quantum mechanics. 

On the other hand in the last decade, starting from 
Czachor's papers [3, HI, EPR correlations in the relativis- 



tic context have been widely discussed [1, 



ic , n, 

_ _.. 25, 26, 

271, 128|, [MEOdMlIlJlS IS , 

43]. Unfortunately the incompleteness of the relativistic 
quantum mechanics formalism (e.q. lack of the covariant 
notion of localization) [52j causes that our understanding 
of relativistic aspects of quantum information theory is 
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far from being satisfactory. Moreover it is unclear which 
spin operator should be used in the relativistic context 
(spin is not a self-contained, irreducible geometrical ob- 
ject in the relativistic quantum mechanics). There are 
several different operators which has the proper nonrel- 
ativistic limit (for the discussion see Sec. IfVj) . Measure- 
ment of quantum spin correlations in EPR experiments 
could help us to decide which spin operator is more ap- 
propriate. In the present paper we use quantum field the- 
ory methods to calculate correlation function for massive 
Dirac particle-anti-particle pair in the EPR type experi- 
ment using particular, in our opinion the most adequate, 
spin operator. Our results can be useful for the discus- 
sion of EPR experiments in which the spin correlation 
function of elementary particles is measured. Such an 
experiment has been performed in seventies [idj but with 
the nonrelativistic particles (protons). 

Entangled pairs of massive particles are usually cre- 
ated in the decay processes of elementary particles (e.g. 
7T° — > e + e~, Z° — > e + e~). Decaying particle has of 
course well defined Poincare-covariant state (e.g. ir° is 
a pseudoscalar particle, Z° four-vector one). Dynamics 
of the decay process is Poincare invariant therefore the 
two-particle entangled state of decay products also posses 
well defined transformation properties with respect to the 
full Poincare group (compare [181]). In our paper we clas- 
sify such states and calculate correlation functions in the 
pseudoscalar and four-vector states which correspond to 
the states of the pair created in the ir° and Z° decay, 
respectively. Our results are valid for any Dirac particle- 
anti-particle pair but we concentrate our discussion on 
the 7T° — > e+e~ and Z° — ► e + e~ decays because particles 
produced in these decays are ultrarelativistic. 

In Sec. |TT] we establish notation and briefly remind ba- 
sic facts concerning free quantum Dirac field. In Sec. IIIII 
we consider two-particle states and classify them accord- 
ing to transformation properties with respect to the full 
Lorentz group. Next section we devote to the discussion 
of the spin operator. Finally, in Sec. |V]we calculate ex- 
plicitly correlation function for the particle-anti-particle 
pair in the pseudoscalar and vector state. The last sec- 
tion contains our concluding remarks. 
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II. THE SETTING 

Let the field operator fulfills the Dirac equation 
(iy*0 M - m)$(x) = 0, (1) 

where 7^ are Dirac matrices (their explicit form used in 
the present paper and related conventions can be found 
in the Appendix [Al) . The field transforms under Lorentz 
transformations according to: 



/7(A)*(a;)J7(A- 1 ) = D(A- 1 )<l(Aa;), 



(2) 



where U(A) belongs to the unitary irreducible represen- 
tation of the Poincare group and D(A) is the bispinor 
representation of the Lorentz group (see Appendix |EJ . 
Field operator has the standard momentum expansion 



*=±k J 



2fc° 



-ikx 



u CXG {k)a a {k) 



x v aa (k)bl(k)], (3) 



where a\{k) (fc)) are creation operators of the particle 
(antiparticle) with four-momentum k and spin compo- 
nent along z-axis equal to a, and k° = %/ m 2 + k. These 
operators fulfill the standard canonical anticommutation 
relations 



{a a ,(k'),aUk)} = 2k°5 3 (k - k')5„ 
{b„,(k'),bl(k)} = 2k°6 3 (k-k')5 a 



(4) 
(5) 



and all other anticommutators vanish. The one-particle 
and antiparticle state with momentum k and spin com- 
ponent a arc defined as 



|fc,<7> o = at(fc)|0>, \k,a) b = bt(k)\0), 



(6) 



respectively. Here |0) denotes Lorentz-invariant vacuum, 
(0|0) = 1, a a (k)\0) = 0. The states © span the carrier 
space of the irreducible unitary representation U(A) of 
the Poincare group 

17(A)!*, a) a/b = T>(R(A, k)) Xa \Ak, A) o/6 , (7) 

where T> £ SU(2) is the matrix spin 1/2 representation of 
the SO(3) group, R(A, k) = L^ALk is the Wigner rota- 
tion and Lk denotes the standard Lorentz boost defined 
by the relations L^k = k, Lr = I, k = (m, 0). As follows 
from Eqs. (HHHJ), states (O are normalized covariantly 



a/b (k,a\k>,<j>) a/b = 2k°6 3 (k~k')5, 



(8) 



Eqs. (|2l6l7p imply standard consistency (Weinberg) con- 
ditions for amplitudes 

v(Ak) = D(A)v(k)V T (R(A, k)), (9) 
u(Ak) = D(A)u(k)V^(R(A, k)), (10) 



where v(k) and u(k) denotes matrices [v^^k)] and 
[u^a^k)], respectively. Explicit form of amplitudes de- 
pends on the chosen representation of gamma matrices. 
Under the choice given in Appendix [A"l amplitudes can be 
written as 



u{k) = iv(fc)<72, 

v(k) = 7 5 v(fc), 



(11) 
(12) 



where, for for the sake of convenience, we have introduced 
the matrix 



v(fc) 



1 



2Jl + £ 



(l+±k^)a 2 
(l+ik^a,)a 2 



(13) 



with k? = {k ,— k), and ai, i = 1,2,3, designate the 
standard Pauli matrices, do = I 53]. 

Action of the charge conjugation C and space inversion 
P on the Dirac field has the form: 



C*(a;)C t = r] C ci> T {x), 



(14) 
(15) 



where ^> — ^7°, \rjp\ = \rjc\ = 1 and C is the charge 
conjugation matrix (|A2[) . 

Operators P and C act on creation operators as follows: 

Pat (fc)pt = Vp aUk*), P&t(fc)P f = -TlpbKk*), (16) 
Ca CT (fc)C t = i lc b a (k), Cftt (fcjd = VcaUk). (17) 

It is convenient to introduce covariant one-particle and 
anti-particle states 

\ a > k )a/b = v a *{k)\k,o) a/b . (18) 

Such states has the following transformation properties 

U(A)\a,k) a/b = D(A- 1 U\(3,k) a/b , (19) 

implied by Eqs. (|7IC6j) . and 

P\u,k) a = rf pl %\p,k«) a , C\a,k) a = v * c \a,k) b , (20) 
P\a,k) b = -VPl ap\P,k*) b , C\a,k) b = r )c \a,k) a . (21) 

III. TWO-PARTICLE STATES 

The main goal of our paper is to calculate spin cor- 
relation function of two Dirac particles. We calculate 
correlations in the states consisting of one particle and 
one anti-particle. The space of such two-particle states 
is spanned by the vectors 



\(k )( 7) b ,(p,X) a ) = bt(k)a{(p)\0). 



(22) 



To analyze transformation properties of the states it is 
convenient to introduce the covariant basis, analogous to 
OH), defined by 

\(a,k),{P,p))=v(k) au v(p)p X \(k,a) b ,(p,\) a ), (23) 
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Notice that for simplicity we have omitted indices a and b 
in the state vector on the left hand-side of (|2"3"|) . Hereafter 
we use the convention that in the two-particle state vec- 
tor left pair of indexes refers to anti-particle, right pair 
to particle, respectively. Using (|C6[) we have 



U(A)\(a,k),((3,p)) 

= D(A- 1 ) aa ,D(A- 1 ) m \(a',Ak),(pi 1 A P )}. 



(24) 



Moreover, we can easily determine the action of the dis- 
crete operations P, C on the states (|2"3"]l 

P\(a,k),(/3,p)) =-i aa ,^ f} ,\(a>,k*),(P',P*)), (25) 
C\(a,k),(p,p))=-\(J3,p) t {a,k)). (26) 

Particle-anti-particle pairs are usually created in the 
state in which the total four-momentum is determined as 
for example in the decay 7r° — > e~e + where the total four- 
momentum of the electron-positron pair is equal to the 
four-momentum of the decaying 7r° . (This decay channel 
has very small but nonzero width [45|.) The most general 
particle-anti-particle state with total four-momentum q 
has the form 



lv> e 



d 3 kd 3 p . 



[V^(fc,p)CrJ \(a,k), (0,p)), (27) 



where ip A (k,p) arc numerical functions, C is given by 
121) and Ta forms a subset of matrices 



(28) 



which transform covariantly under Lorentz transforma- 
tions and form a basis of the Clifford algebra generated 
by Dirac gamma matrices. Matrices T a transform with 
respect to the index A according to certain representation 
D(A) of the Lorentz group 



D(A)r A D(A- 1 ) -B AB (A^ 1 )r B . 



(29) 



Note that we have inserted matrix C in Eq. (|2"Tj) because 
from {Ml we have D T (A)C = CD _1 (A). 

Identification of the singlet or vector state is based 
on the transformation properties of the state (|27p under 
Lorentz transformations and parity P. Using (|24j) we have 



A 



^^' A (k,p)Cr A \(a,k),(fi,p)), (30) 



a/3 



where from Eq. (|29[) we obtain 



M (fc,p) = P^(A- 1 )^ b (A- 1 A : , A-'p). (31) 



Under parity the elements of the Clifford algebra trans- 
form as follows 



1 Ta7° 



ABt B, 



(32) 



where P represents parity in the carrier space of the rep- 
resentation D; namely, the state transforming accord- 
ing to {22!) and (02]) will be called: scalar if D(A) = I, 
P = I, {T A } = {I}; pseudoscalar if D(A) = I, P = -I, 
{T A } = { 7 5 }; four-vector if D(A) = A, P = rj = 
diag(l, — 1, — 1, — 1) , {r^} = {7^} and pseudo-four- 
vector if D(A) = A, P = -r), {T A } = {7 M 7 5 }- 

In this paper we restrict ourselves only to the pseu- 
doscalar and four-vector case since such states can be 
identified with real particles decaying into e + e~ pair 
(compare for ex amp le ir° and Z° and their corresponding 
decay channels |45j). 

a. Pseudoscalar The general pseudoscalar state 
reads 

\< = jw^ {q ~ k ~ pMKp) 

x (C 7 5 ) a/3 |(a,fc),GM> 

x (v T (fc)C 7 5 v(p)) \(k,a),(p,X)), (33) 

where, according to Eq. (f3"Tj) the density function trans- 
forms under Lorentz transformations as follows: 



<p'(Ak,Ap) = <p(k,p). 



(34) 



b. Four-vector The most general four-vector state 
has the following form 

x (v r (fc)C 7 „v(p)) \(k,a),{p,X)). (35) 

\ / crA 

One can easily check that 

(Pn + k tl )v T {k)C^v( P ) = 0. (36) 



Therefore, in the integral (|35|) the nonzero contribution 
has only such Lp^(k,p) : which fulfills the transversality 
condition 



q^(k,p) = (p M + k„)^(k,p) = 0. 



(37) 



Notice that in the decay cases the functions ip(k,p), 
(p^(k,p) are related to the dynamics of the decay. 



IV. RELATIVISTIC SPIN OPERATOR 

To calculate correlation function we have to introduce 
the spin operator for a relativistic massive particle. In the 
discussion of the relativistic Einstein-Podolsky-Rosen 
experiment various authors use different spin operators 
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0, 0, HI SI H HE H||. However it seems [H, [13, HI 
that the best candidate for the relativistic spin operator 



S = — W 
m \ 



(38) 



where is the Pauli-Lubanski four-vector 

W" = \e^ aX Pj aX . (39) 

Here P v is a four-momentum operator and J„\ de- 
notes generators of the Lorentz group, i.e., U(A) — 
exp (iuj^ J^ v ). One can show [46[ that operator ([38]) is 
the only operator which is a linear function of W^ L and 
fulfills the relations 



[J\S j ] =ie ijk S k , 
[P",#]=0, 



(40a) 
(40b) 
(40c) 



and is a pseudovector, i.e. PS l P = S\ Here J 1 



In the representation (|A1|) of gamma matrices we have 



W°\(a,k))=W° f3 \((3,k)), 
W\(a,k))=W a0 \((3,k)), 



where 



1 fk-a- 

2 I kcr 



W = — 



fc°<X 



/ 
/ 



i(k x <t) ( _° 7 



Therefore for the spin operator we easily find 
1 



S|fc,cr) 
S|a, k) 



-axa\k,X), 



v(fc)^-v(*) 



(41) 
(42) 



(43) 
(44) 

(45) 
(46) 



a/3 



Real detectors register only particles which momenta be- 
long to some definite region in momentum space. There- 
fore, taking into account Eqs. (|40cl45[) . in one-particle 
subspace of the Fock space, the operator measuring spin 



of the particle with four-momentum from the region A 
of the momentum space has the form 



d 3 k 

2k° 



Uk)^a x (k). 



(47) 



This operator gives zero when acting on the anti-particle 
state or state of the particle with four-momentum out- 
side the region A. Analogous operator measuring spin of 
the anti-particle with four-momentum belonging to the 
region B of the momentum space can be written as 



fk 

2fc° 



bt(k)^b x (k). 



We have 



Ml _ 



0. 



(48) 



(49) 



V. CORRELATION FUNCTIONS 



In this section we calculate correlation function in the 
EPR-Bohm type experiment. In such an experiment we 
have two distant observers, say Alice and Bob. We as- 
sume that both observers are at rest with respect to the 
same inertial reference frame. Particle-anti-particle pair 
is produced in an entangled state \(p), the particle is regis- 
tered by Alice while the antiparticle by Bob. Alice mea- 
sures spin component of the particle along direction a, 
Bob spin component of the anti-particle along direction 
b. Therefore Alice uses the observable a-S^ and Bob the 
observable h-S b B (see Eqs. (1471481) V So, by virtue of ([4^]) . 
the normalized correlation function in the state \ip) has 
the following form: 



C^ B (a,b)=4 



(<p\ (a-S^) (hty) \<p) 



(<P\<P) 



(50) 



We calculate correlation function in two important cases 
when EPR pair is produced in the pseudo-scalar or four- 
vector state. 



A. Pseudoscalar state 

The pseudoscalar state is given by Eq. (|3"3"1) . Therefore 
from (l50l we find 
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C£ B ( a ,b) = 



d 3 k d 3 p 

2fc° 2^ 



lxA{p)XB{k) {5% - k-p)f |^(fc,p)| 2 Tr b. < T(v T (fc)C 7 5 vb))a.<T T (v T (fc)C7 5 v(p)) t 



(51) 



where xa{p) i Xb(&) are characteristic functions of the regions A and B in the corresponding momentum spaces. 
Using Eqs. (fl~3l IA2jl and (|Alj) we arrive after little algebra at 



Tr 



b-(T v' (fc)C 7 5 v(p) a-er J v J (fc)C 7 5 v(p) 



— (a-b(m 2 + fcp) - (k x p) 



a x b 



(b-p)(a x k) - (ak)(b x p) 



(m + fc°)(m + p°) 



}, (52) 



Tr 



(v T (fc)C 7 5 v(p)) (v T (fc)C 7 5 v(p) 



m + kp 



(53) 



Therefore in the simple situation when momenta of both particles in the state \ f) q are sharp (i.e. characteristic 
functions in Eq. (|5ip are replaced by delta functions) we obtain the following correlation function: 



C£(a,b)=-a-b- 



k x p 



m 2 + kp 



a x b 



(b-p)(a x k) - (ak)(b x p) 

(m + k°)(m +p°) 



(54) 



In this special case of the sharp momenta we obtained the same correlation function (|54[) in our previous paper [l2| 
where we discussed the Lorentz-covariant spin density matrix in the framework of the relativistic quantum mechanics. 
Notice that in the case when the laboratory frame (observers) coincides with the center of mass frame, or even when 
p||k, from Eq. ([54")) we get the same correlation function as for the singlet in the nonrelativistic case 



C p c 5 MF (a,b) = -ab = C p k 5 IIP(a,b). 



(55) 



It is interesting that correlation function calculated by Czachor in [5j, in the center of mass frame still depends on 
momentum and has the following form 



a-b — 1 a^-bj_ 



Vl + ^[(n-a)2-iyi + /?2[(n.b)2-l]' 



(56) 



where n = k/|k|, /? = |k|/fc°, a^ = a — (n-a)n and b_L = b (n-b)n. 

Formulas (|55f and (|56[) are different because Czachor uses different spin operator. Therefore the experimental 
measurement of the correlation function in the center of mass frame could show which spin operator is more adequate 
in the relativistic quantum mechanics. One of the possible sources of the electron-positron pairs is the n° decay into 



the channel ir° 



Since m n o = 134.98 MeV/c 2 , 



0.51 MeV/c [45J, electrons and positrons produced in 



this decay in the center of mass frame are ultrarelativistic. Therefore let us find the limit /3 
and ([56]>. We get 



1 of the formulas (155 



,-,CMF 
°Czachor 



C p ? F (a,b) 
( a > b ) \a-+i = 



= ab, 



(n-a)(n-b) 
l(n-a)(n.b)| 



±1, 



(57) 
(58) 



respectively. We point out the discontinuity in the Czachor's correlation function in the ultrarelativistic limit (|58[) . It 
is also interesting to notice, that the function 



AC = C c ^ho r (a,b)-C p - s ^(a,b) 

can take quite large values for (3 close to 1. We shown this function in the Fig. Q] for the value (3 
corresponds to e + e~ created in the ir° decay at rest. 



(59) 

0.999 which 



* „ 

FIG. 1: Function AC (Eq. (|59[l ) in the parametrization n = (0, 0, 1), a = (sin 9 a , 0, cos 8 a ), b = (cos (pi sin 9b, sin (pi sin Oi,, cos 8t). 
We have plotted the graph of the function AC for two arbitrarily chosen values of tpi, and for ft = 0.999, which is the value for 
electron-positron pair created in the 7r° decay. 



B. Vector state 



General vector state is given by Eq. (|35j) . So in this case from (|50|) we get 

C*ve? ( a j k) = 

' /k ^JuOOxbW (<5 4 (<z - t-pf^^pK^P^ [b.<r(v T (A : )C7' J v(p))a. < T T (v T (fc)C7 i 'v(p) 



2fc° 2p° 



^^XA(p)XB(k) (S^q k p))\,(k,p)rt(k,p) Tr [(v T (fcXVv(p)) (v T (fc)C 7 i/ v(p)) t 



(60) 



where, as before, Xa(p) 1 XB{k) are characteristic functions of the regions A and B in momentum space. If the 
observers frame coincides with the center of mass reference frame (in which p = k n = (k , — k)) then by means of the 
transversality condition (|37p , we have 



ha ( v T (fc)C 7 ''v(fc 7r ) ) a<r T ( v T (fc)C 7 l/ v(fc 7r ) 



-U (a-b) ffcO", ,a _ | k pi _ ^ [(^(b^) + (a^b-*,*)] - 2 ( a ' k )( b - k )l k ^ 2 

771 L J 



fc° 



(m + fc ) 2 

(a-k) ((b-^^-k) + (b- V *)( V -k)) + (b-k) ((a-^^-k) + (a. V *)(yk)) [>. (61) 



^(M"V:(Wrr 



(v T (fc)C 7 ^v(fc")) (v T (A;)C 7 I/ v(A; ,r )) f 



7/7 " 



fc°>| 2 -|k-vl 



(62) 



Therefore, if we assume that momenta of the particles in the state \fY^ are sharp, the correlation function in the 
center of mass frame are given by (antiparticle has momentum k, particle fc 7r ) 



k° 



Crr(a,b) = a b 2 

fc° 2 |^| 2 - |k-^| 2 



[(aV)(b-<p) + (a-^Cb.^*)] - 



2(a-k)(b-k)|k-<p| 



{m + k°) 2 (k° 2 \<p\ 2 - |k-<p| 2 ) 
k° [(a-k) ((b-y)(k-y») + (b-y*)(k-y)) + (b-k) ((a-y)(k-y») + (a-y*)(k-g))] 
{m + k°)(k° 2 \ip\ 2 -\k-<p\ 2 ) 



(63) 
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Notice that in this case, in opposition to the one considered previously, in the center of mass frame the correlation 
function depends explicitly on the momentum. Moreover, if k _L <p from (|63[) we get the same result as for nonrela- 

tivistic triplet state (|D5[) . Also in the nonrelativistic limit (k° 2 3> |k| 2 ) from (|63[) we get the correlation function for 
nonrelativistic triplet state (|D5|1 . 

When we consider electron-positron pairs produced in the Z° decay two remarks are in order. First, in some 
experiments it is possible to produce polarized Z° (SLAC, Stanford [471] ). so in such experiments the vector (p is 
determined. Second, electron and positrons produced in the channel Z° — > e + e~ are highly ultrarelativistic. In the 
limit (3 — > 1 the formula (f6"3")) takes the form 

a°ec F (a,b)|^ 1 = a-b - | | 2 {( a -y*)( b -y) + (*¥>)(*>■¥>*) + 2(a.n)(b.n)|n.^| 2 

- (a-n) [(b-<p)(n-cp*) + (b-tp*)(n-tp)] - (b-n) [( a -ip)(n-<p*) + (a-cp*)(n-<p)] }, (64) 

where n = k/|k|. Eq. (|64p takes very simple form in the configuration a _L ip, b _L cp: 

rMf , N , 2(a-n)(b-n)|n-<p| 2 , . 

C™ F (a ± v ,b± <p)\^ x = a-b - , (65) 

It is interesting to notice, that the function 

AC vec = C™ F (a l Vl bl <p)L^ - C nonrelativ = - 2i * n l (h * )ln :f , (66) 

where C non reiativ denotes correlation function for the nonrelativistic triplet state (ID5|i . can take arbitrary value from 
the interval [—2,2]. We shown this function in the Fig.O 

I 




FIG. 2: Function AC (Eq. (|66jl ) in the parametrization 
(p = (0,0,1), n = (sin#, 0, cosfl), a = (cos(^ a , 8mtp a , 0), 
b = (cos ifb, sin ipt,, 0). In this parametrization AC = 
— 2 cos if a cos cos 2 We have plotted the graph for the 
value 6 such that cos 2 6 — 1. We see that AC is of the same 
order of magnitude as the correlation function. 



VI. CONCLUSIONS 

We have discussed correlations of two relativistic mas- 
sive particles in the EPR-type experiments in the context 
of quantum field theory formalism. Choosing the most 
appropriate spin operator (|38p , we have calculated corre- 



lation function for the particle-anti-particle pair in the 
pseudoscalar and vector states. We have found general 
formulas and discussed them in some details for the sharp 
momentum states in the specific configurations. Rela- 
tivistic correlations in the vector state have been never 
discussed before. It should also be noted that all these 
functions posses a proper nonrelativistic limit. 

It is interesting that for the pseudoscalar (singlet) state 
the correlation function in the center of mass frame is 
the same as in the nonrelativistic case. Therefore, in 
opposition to Czachor's results the degree of violation 
of Bell inequality for the pseudoscalar state in CMF is 
independent of the particle momentum. The correlation 
function for the vector state depends on the configuration 
and for same configurations gives the same result as for 
the nonrelativistic triplet state. 

It seems that EPR-type experiment with relativistic 
elementary particles are feasible with the present tech- 
nology |48j . 
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APPENDIX A: DIRAC MATRICES 

In this paper we use the following conventions. Dirac 
matrices fulfill the relation 7^7" + 7 I/ 7 M = 2<7 MI/ where 
the Minkowski metric tensor — diag(l, —1,-1,-1); 
moreover we adopt the convention g 0123 = 1 . We use the 
following explicit representation of gamma matrices: 

M?o)' -UT). "Mo- ,)' 

A 1 1 

where <x = (01, a%, 03) and <7, are standard Pauli matri- 
ces. The charge conjugation matrix has the form 

C = ^ 2 1 Q = l (%> _°J, (A2) 

so 

APPENDIX B: BISPINOR REPRESENTATION 

By the bispinor representation of the Lorentz group we 
mean the representation D^^) ® £)(°>i). Explicitly, if 
A e SL(2, C) and A(A) is an image of A in the canonical 
homomorphism of the SL(2, C) group onto the Lorentz 
group, we take the chiral form of ®D(°'i\ namely 

D(A(A)) = (jj { J r ^j . (Bl) 

The canonical homomorphism between the group 
SL(2, C) (universal covering of the proper ortochronous 
Lorentz group Ll_) and the Lorentz group L, ~ 
50(l,3)o [49] is defined as follows: With every four- 
vector k^ we associate a two-dimensional hermitian ma- 
trix k^a^, where (Tj, i = 1,2,3, are the standard Pauli 
matrices and cto = /. In the space of two-dimensional 
hermitian matrices the Lorentz group action is given by 
k'^CTfj, = A(k fJ 'a t _ t )A^ , where A denotes the element of the 
5L(2, C) group corresponding to the Lorentz transfor- 
mation A(A) which converts the four- vector k to k' (i.e., 
= A» v k v ). 

APPENDIX C: USEFUL FORMULAS 

The following relations hold: 



(k-f) v(k) = mv(fc), 


(CI) 


v(fc)v(fe) = 7, 


(C2) 


7 °v(fc) = v(fc-), 


(C3) 






v(fc) 7 ^v(fc) = —J, 
m 


(C4) 


v ( fe M*) = 37^7 + wl), 


(C5) 



and 

v(Afc) = D(A)v(fc)P T (i?(A, fc)), (C6) 
v(Ak)a 2 = D(A) v(/fc)(T 2 V\R(A, k)). (C7) 

Notice that when operator ^4 acts on standard basis vec- 
tors in the following way 

A\k,a)=A aa ,\k,a'), (C8) 

then its action on the covariant basis (fT5)) is of the form 

A\a,k) = {y{k)Ay{k)) aoi ,\a',k). (C9) 

Analogously for two-particle states the relation 

(i ® B)\(k, a) b , (p, A) Q ) = A CTCT ,5 AA , |(fc, </)b, (p, 

(CIO) 

implies 

(1® £)](«, fc), CM) = (v(fc)iv(fc)) 

(v(p)j}v(p))^|(a',fc),(/3',p)). (Cll) 

APPENDIX D: CORRELATIONS IN 
NONRELATIVISTIC TRIPLET STATE 

We recall in this appendix the correlation function 
in nonrelativistic triplet state (see e.g. [50]). Let |cr), 
a = ±5, denotes eigenvector of the spin component along 
z axis corresponding to the value of spin z-component 
equal to a . The general triplet state has the form 

W)= ]T p„Ak> ® |a>, (di) 

cr,A=±l/2 

with the symmetry condition 

<AxA = ¥>Acr- (D2) 

It is convenient to parametrize matrix ip = [p a \] in the 
following way: 

cp = —(<p.o)a2. (D3) 

So finally the normalized correlation function in the 
triplet state is given by 

C(a,b) = (<p|a-er® b-cr\ip) (D4) 
= a-b - -L^Xb^*) + (b- V >)(a- V >*ip5) 
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